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NOTATION 


Wave  amplitude 

Added  mass  coefficient  in  the  ith  mode  due  to  motion  in  the  jth  mode 
Waterplane  area  at  load  waterline 
Sectional  heave  added  mass 

Damping  coefficient  in  the  itli  mode  due  to  the  motion  in  the  jth 
mode 

Vertical  force  rate  due  to  heave  velocity 
Vertical  force  rate  due  to  pitch  angular  velocity 
Pitch  moment  rate  due  to  pitch  angular  velocity 
Pitch  moment  rate  due  to  heave  velocity 

Vertical  distance  between  the  center  of  buoyancy  and  the  center  of 
gravity  of  ship 

Sectional  heave  damping  coefficient 

Restoring  coefficient  in  the  ith  mode  due  to  the  motion  of  the  jth 
mode 

Lift-curve  slope  of  the  ith  fin 
Chord  of  the  ith  fin 

Wave  exciting  coefficients  us  defined  below  Equations  (14) 

Longitudinal  metacentric  height 

Gravitational  acceleration 

Mass  moment  of  inertia  of  ship  about  the  y  axis 

Overall  ship  length 

Length  between  Stations  0  and  20 

The  x  coordinate  of  the  1 /4-chord  point  of  the  ith  fin 

Mass  of  ship 

Coefficient  of  pitch  moment  due  to  pitch  angular  velocity 
Coefficient  of  pitch  moment  due  to  heave  velocity 
Moment  of  waterplane  area  about  the  y  axis 
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Right-handed  Cartesian  coordinate  system;  see  definition  on  page  2 
Span  of  the  ith  fin 

Ship  speed 

Coefficient  of  heave  force  due  to  pitch  angular  velocity 
Coefficient  of  heave  force  due  to  heave  velocity 

Displacement  of  ship 
Argument  of  liw  -  X, )  (iw  -  X2 ) 

Argument  of  <iw  -  X} )  <iu>  -  X4 ) 

Damping  ratio,  i  =  3  for  heave  and  i  *  5  for  pitch 
(=XnR  +  i  Xnl)  nth  stability  root 
Heave  displacement,  positive  upward 
Pitch  angular  displacement,  positive  bow  down 
Mass  density  of  water 

Wave  encounter  frequency  in  radian  per  unit  time 
Undamped  natural  frequency  for  heave  mode 

Undamped  natural  frequency  for  pitch  mode 
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INTRODUCTION 

r: :;;“™  rr  rr  —  --  -  -  - 

Tht  Period,  together  with  smaller  "*  '****' 

merged  hulls,  provides  SWATH  chine  ^ „  ,  xtitation  lorec  due  to  the  sub- 

smallness  of  the  waterplane  area,  hMJJ.tr^^'^.h-.,,,0,fc”,e 

cruises  with  speeds  higher  than  a  certain  limit  The  cause  of  Mi  '  "l  h'"  “  ihip 

stems  from  the  so-called  Monk  moment,  which  is  proportional  mai"ly 

end  provides  a  destabilising  pitch  moment.  Fo,  the  SWATH  4  1  M  T™  "  “*  ^ 

Hus  instability  would  occur  at  appioximalely  ;7  knots  ‘  W“  PrC'"t,"J 

,  r  rr fins 

coupled  pitch  and  heave  stability  conations  ,i„  e  ,  '  ^  CTm'"ed  ™  «*  W  of 

One  fin  was  installed  on  the  inboard  side  of  each  huTof 'sWATO  r™1  '°  '* 

position.  Free-running  model  experiments  in  waves  were  carried  ollT'T^  °**  1 
the  stabilizing  fins.  The  results  of  mmini  ho  1  w,t  1  a,ul  w',llont 

indeed  occur  when  the  full-scale  speed  exceerrCknor^nTth'1'31  dkl 

:  :• — -  «*- 

Z'Z’z:  ri: ■ :  r  "•  ~ »™  -■  “  r.. 

0.833  to  I  while  SWATH  4B  ,,''  *  °f  WATO  4  *»  ■  -ducecl-scalc  ratio  of 

I.  wnne  bWA  fH  4B  has  the  same  displacement  as  su/atu  a  a  u 

length,  a  smaller  hull  diameter,  a  greater  strut  length  a  H  4A  ^  °VeraH 

4B  has  a  longitudinal  metacentric  height  (C  M  )  •  b  "  Sma,,er  S‘n,t  beam-  Thc  SWA™ 

temr.c  might  (C,MC)  about  twice  as  large  as  SWATH  4A. 
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A  SWATH  >hip  with  two  separated  struts  at  fore  and  aft  can  provide  a  large  GMg  It 
is  of  interest  to  learn  the  effects  of  increased  GMg  on  the  vertical  plane  stability  as  well  as 
on  the  vertical  motion  in  waves  in  the  vicinity  of  resonant  frequencies.  The  reason  for  choos¬ 
ing  SWATH  4B  is  intended  to  pursue  the  aforesaid  interest.  Figure  I  shows  the  profiles  of 
SWATH  4 A  and  SWATH  4B.  and  Table  I  shows  the  principal  characteristics. 

The  analysis  for  determining  adequate  fin  sizes  is  based  on  the  coupled  heave  and  pitch 
equations  for  stability  and  will  be  described  in  detail  in  the  following  sections.  The  size  of 
the  fins  is  determined  by  the  following  three  conditions.  The  fins  must  provide: 

1 .  pitch  stability  up  to  40  knots  in  calm  water. 

2.  increased  damping  for  the  heave  and  pitch  modes,  and 

3.  reasonable  natural  periods  for  heave  and  pitch. 

DESCRIPTION  OF  ANALYSIS 

Assume  that  a  SWATH  ship  translating  with  a  forward  speed  U  is  momentarily  disturbed 
so  as  to  induce  a  coupled  heave  and  pitch  motion.  Our  objective  is  to  study  whether  the 
ship  at  a  given  speed  has  sufficient  stability  to  restore  itself  to  its  original  equilibrium  in  a 
reasonable  time.  If  the  ship  is  found  unstable,  a  stabilizing  fin  or  fins  will  be  introduced  in 
the  analysis. 

The  coupled  heave  and  pitch  equations  of  motion  without  excitation  can  be  expressed 
in  the  form 

(M  ♦  A,,*,  *  B),t,  ♦('„(,  -  Ajjtj  +  Bjjlj  +  Cjjtj  =0  (II 

a,,!,  *  b5)  i,  ♦  i-„{,  ♦  <i,  *  a„ )t5  ♦  Ms  +c„is  -o  ■:> 

Equations  (I )  and  (2)  are  formulated  with  respect  to  the  body  coordinate  system  at  the 
equilibrium  position  of  the  ship,  i.e.,  no  motion  except  a  steady  forward  translation.  The 
origin  of  the  coordinate  system  0  is  located  on  the  calm  water  free  surface  and.  together  with 
the  center  of  gravity,  lies  in  the  longitudinal  plane  of  symmetry  of  the  two  hulls.  The  axes. 
Ox.  Oy,  and  Oz  are  respectively  directed  toward  bow.  port,  and  vertically  upward.  $3(t)  is 
the  heave  displacement,  which  is  positive  upward,  and  $5(t)  is  the  pitch  angular  displacement 
about  the  y-axis.  which  is  positive  bow  down.  M  is  the  mass  of  the  ship;  I5  is  the  mass 
moment  of  inertia  about  the  y-axis;  A^,  Bjr  and  Ctj  for  ij  *  3  and  5  are.  respectively,  the 
added  inertia,  damping,  and  restoring  coefficients. 
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TABLK  I  HULL  ('HARACTHRISTICS  OF  SWATH  4A  AND  SWATH  4B 


Description 

Unit 

SWATH  4 A 

SWATH  48 

Displacement 

Long  tons 

7379 

7379 

Displacement  Mam  Hull 

Long  tons 

IBS? 

1829 

Length  Overall 

Feet 

739  7 

775  9 

Length  of  Strut 

Feet 

188  9 

779  54 

Minimum  Strut  Thickness 

Feet 

6  67 

5  147 

Maximum  Hull  Diameter 

Feet 

ISO 

13  7 

Draft 

Feet 

25  0 

7?  63 

Waterplane  Area 

Feet’ 

1893 

1975 

Water  plane  Area 

Moment  of  Inertia 

Feet4 

3  645  *  10*’ 

5  997  a  106 

Center  of  Buoyancy 
from  Nose  of  Hull 

Feet 

1156 

133  92 

Center  of  Flotation 

From  Nose  of  Hull 

Feet 

117  7 

135  93 

BG 

Feet 

15  19 

14  06 

Longitudinal  GM 

Feet 

77  7 

58  4 

Radius  of  Gyration 
for  Pitch 

Feet  | 

5?  1 

61  8 

Figure  1  -  Profiles  of  SWATH  4A  and  SWATH  4B 

(Number i  are  in  feel,  and  drawtnis  are  nol  to  scale) 
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It  will  be  assumed  that  the  hydrodynamic  coefficients  AtJ  and  B(j  are  made  up  of  three 
contributing  factors 

I  irrotationa!  (low  with  the  free-surface  boundary  disturbed  by  the  body  without  the 
fin, 

2.  real  fluid  effect  on  the  bare  hull,  and 

3.  fin  effect. 

Under  this  assumption  the  hydrodynamic  and  restoring  coefficients  are  given  by 

Ajj  =J  a33<xkix  (3,I) 

L 

,  , 

BJ3  *  J  b33,x)dx  +  UB33  7PUC'S'C“*'  <3‘2) 

i*  1 


C33  =  P«Awp  +  C33* 


*35  = 


*  j  Xa33 


(xKix 


C„  -  -  «M„P  ♦  UB„  +  C% 
Aj,  j  sijj 


,(x)dx 


Cj,  =  -  pgMwp  +  CJ3 

ASJ  =  J  x2a33 


.(x)dx 


(3.3) 

(3.4) 

£-»  1 

^33  +  B?5l-2T^UcMCLai 

1*  1 

(3.5) 

-(F) 

-35 

(3.6) 

(3.7) 

1 

53  T^^Cisi®i^Lai 

i  1 1 

(3.8) 

(3.9 

(3.10) 

4 


(3.11) 


0  1^ 

B5S  =J  x^jjUWx  +  UB*,  vtfciai 

i“  I 

C5i  =  A  GMg  +  ^pU  W  W ‘iV 

I*  I 

where  I  =  integration  over  the  ship  length 

a3J(x)  =  sectional  added  mass  of  cross  section  at  x 
h,,(x)  =  sectional  damping  due  to  wavemaking  of  cross  section  at  x 
B3  3  =  ver,'cal  force  rate  due  to  heave  velocity* 

N  =  total  number  of  fins 

P  =  mass  density  of  water 

Cj  =  chord  of  the  ith  fin 

s,  =  span  of  the  ith  fin 

CLai  =  lift-curve  slope  of  the  ith  fin 

g  =  gravitational  acceleration 
Awp  =  waterplane  area 

B*j  =  vertical  force  rate  due  to  pitch  angular  velocity* 

=  x  coordinate  of  the  ith  fin 
MwP  =  waterplane  area  moment  about  the  y  axis 
Bjj  -  pitch  moment  rate  due  to  heave  velocity 
Bss  =  P'tt-h  moment  rate  due  to  pitch  angular  velocity 
A  =  displacement  of  ship 
GMg  =  longitudinal  metacentric  height 

( M  I 

C55  =  Munk  moment  =  B£3 

Tlie  coefficient  C  ■■  for  i,  j  =  3  and  5  are  to  account  for  the  effects  of  sinkage  and  trim 
resulting  from  forward  speed  of  the  ship  in  the  free  surface  on  the  heave  force  and  pitch 


•Contribution  of  bare  hull  only  and  refers  to  nonwavemaking  part. 
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moment  derivatives.  The  effect  of  these  coefficients  on  the  stability  analysis  were  assumed 
to  be  negligible.  The  importance  of  these  terms  increases  as  the  draft  decreases. 

The  coefficients  with  an  asterisk  represent  the  contributions  of  the  hull  without  fins  to 
the  force  and  moment  not  attributable  to  wavemaking.  Since  these  quantities  are  strongly 
influenced  by  viscous  effects,  it  was  necessary  to  estimate  them  on  the  basis  of  available 
experimental  data  from  experiments  with  submarine  models.  The  viscous  effects  thus  obtained 
ignore  the  presence  of  the  struts  and  also  the  hydrodynamic  interference  effect  from  the  other 
hull.  If  we  use  the  notation  employed  in  the  equations  of  motion  for  submarines,1  which  are 
often  called  submarine  derivatives,  then  we  may  write 

B?)  ■  4  PLJz; 

■  -  T 

(3.13) 

H?3  =  ~  pLX 

Bjs  ■  -  4  plX 


where  L  is  the  overall  length  of  the  sliip. 

Although  it  is  known  that  several  of  the  coefficients  of  Equations  (1)  and  (2)  are  fre¬ 
quency  dependent,  they  have  been  treated  as  constants,  and  the  values  selected  were  those 
associated  with  frequencies  close  to  the  natural  heave  frequency.  It  was  felt  that  this  was  a 
reasonable  approximation  for  the  following  reasons. 

1 .  Since  SWATH  vehicles  are  more  deeply  submerged  than  surface  craft,  the  frequency 
dependence  of  the  coefficients  are  significantly  smaller;  especially  the  frequency  range  about 
the  natural  heave  frequency. 

2.  For  SWATH  vehicles  of  this  study,  the  oscillatory  response  to  disturbances  is  domi¬ 
nated  by  motions  at  a  single  frequency,  viz.,  the  heave  natural  frequency.  The  added  mass 
terms  vary  only  mildly  with  frequency  here.  Furthermore,  these  terms  represent  only  about 
40  percent  of  the  virtual  mass. 


*Gertler,  M.  and  G.  R.  Hagen, 


"Standard  Equations  fot  Motion  for  Submarine  Simulation,"  NSRDC  Report  2510  (1967). 
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3.  The  purpose  of  the  study  was  to  investigate  the  relative  merits  of  various  stabilizing 
fin  sizes  in  providing  stability  and  adequate  damping  of  the  motions.  The  effects  of  these 
fins  were  generally  much  larger  than  the  effects  due  to  the  terms  that  would  tend  to  vary 
most  with  frequency,  viz.,  the  damping  terms  due  to  wavemaking. 

4.  Comparisons  of  theory  with  data  for  most  of  these  frequency-dependent  damping 
coefficients  often  show  different*. ..  of  the  same  order  of  magnitude  as  the  variations  with 
frequency. 

Equations  (I)  and  (2)  are  a  pair  of  linear,  homogeneous  differential  equations  with 
constant  coefficients.  Solutions  of  such  equations  are  well  known2  and  are  given  by  the  real 
part  of 


a-‘K' 

n  =  I 

(4) 

N 

n  =  I 

where  an  and  bn  are  arbitrary  constants  that  depend  on  the  initial  conditions,  and  the  X 
are  the  roots  of  the  characteristic  equations.  The  characteristic  equation  may  be  obtained  by 
substituting  any  pair  of  the  solutions  in  Equation  (4)  into  Equations  (1 )  and  (2).  Thus,  if 
we  substitute  the  solutions 


o'1"' 

(5) 

V  ,  X  i 
$5<t)=bn  e  " 

into  Equations  ( I )  and  (2).  we  obtain 

|(M  +  A33)X2  +  B33X  +  C33J  a„  +  (A3jX2  +  B3JX  +  C35)bn  =0 
(A53X2  +  BJ3X  +  C53  )an  +  {tls  +  A55)X2  +  B5JX  +  C5S  }  b„  =  0 


2 Hildebrand.  I  B.. 


"Advanced  Calculus  for  Engineers,”  Prentice  Hall,  Inc.,  Englewood  Cliffs,  N.  J.  (1949). 
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To  obtain  solution*  for  a„  and  b„  from  simultaneous  homogeneous  equations  as 
determinant  of  the  equations  should  vanish  That  is. 


shown,  the 


aX4  +  bXJ  +  cX2  +  dX  +  e  =  0 


(6) 


or,  on  dividing  through  by  a 


X4  +  b’X3  +  c'X2  +  d'X  +  e'  =  0 


where 


a  =  (M  +  A3j)  (l5  +  A55)  -  AjjAjj 
b  =  IM  s-  A,, >BjS  -  AsjBjS  ♦  B„(ls  *  Asi)  -  B„AJS 


(7.1) 

(7.2) 


c  =  (M  + 


a33  yc 


55 


A53C35 


B33B55 


B53  B35 


+  (1 


55 


A55)C 


33 


^35^53 


(7.3) 


d  =  B33C55  “  B53C35  +C33B55  "  C35B53 
e  =  C33C55  "  C35C53 


(7.4) 

(7.5) 


„  is  clear  from  Equation  (6)  that  there  are  four  roots.  These  roots  can  be  obtained 
written  as 


(8.1) 


(8.2) 


3 Dickson,  L.  E.,  “Theory  of  Equations."  John  Wiley  and  Sons,  Inc..  New  York  (1947). 
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1/2 


where  m  = 


c  +  t 


q  =  -  b'V  +  4c'e'~  d'2  +  -y  <b’d'  -  4e')  -  ^  e'3 


As  can  he  seen  from  Equations  (8).  the  value  of  Xn  may  be  real  or  complex,  with  the 
complex  value  appearing  in  conjugate  pairs.  The  motion  corresponding  to  each  real  \  or 
each  complex  pair  is  called  a  normal  mode.  There  are  four  possible  kinds  of  normal  modes 
of  motion  according  to  whether  Xn  is  real  or  complex*  and  whether  it  has  a  positive  or  nega¬ 
tive  real  part.  As  can  be  inferred  from  Equation  <4a)  in  the  footnote  the  motion  is 

1.  divergent,  if  Xn  is  real  and  positive; 

2.  divergent  oscillatory,  if  Xn  is  complex  and  has  a  positive  real  part; 

3.  convergent,  if  Xn  is  real  and  negative;  and 

4.  convergent  (or  damped )  oscillatory,  if  Xn  is  complex  and  has  a  negative  real  part. 
Modes  I  and  2  are  unstable  modes  since,  as  can  be  seen,  the  amplitude  of  the  motion  grows 
with  the  time  as  long  as  the  system  remains  linear.  Modes  3  and  4  are  stable  modes  since 
amplitude  of  the  motion  decays  with  time  to  zero. 

The  necessary  and  sufficient  condition  for  stability,  i.e..  that  all  of  the  roots  of  Equa¬ 
tion  (6)  have  negative  real  parts  is  given  by  Routh's  stability  criteria4 


•If,  for  example,  the  root*  for  one  of  the  modes  were  given  by  Xj  ■  X#R  +  iX^,  and  Xj  ■  XgR  -  iX^,.  then  by  substi- 
tuting  these  into  Equation  (4)  with  N  •  2,  and  taking  the  real  part,  we  would  obtain  the  following  typical  form  of  solution 
for  this  mode 


X_t 

|j(t)  *  c  **  <c(  cos  X^t  ♦  da  sin  X^,t» 


(4a) 


where  X^j  is  the  circular  oscillation  frequency,  and  and  dg  are  real  constants. 

4Routh,  E.  J„  “Dynamics  of  a  System  of  Rigid  Bodies,"  6th  Edition,  The  Macmillan  Company,  London  (1905). 
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a,  b.  d.  e  >  0 
bed  -  ad*  -  b2e  >  0 


(9) 


Numerical  values  of  primary  interest,  which  are  related  to  Xn .  are  the  period  of  oscilla¬ 
tion  Tn.  the  damping  ratio  fn.  and  the  time  to  halve  or  double  the  amplitude  of  the  disturbed 
motion. 

If  the  modes  are  oscillatory,  the  natural  period  are  given  by 

T,  *  2w/lX,,l  (10.1) 


T3  =  2»/|X„l 


(10.2) 


if  we  write  Xn  as  follows 


=  \iR 


iX 


n! 


According  to  Equation  (4).  the  disturbed  motion,  in  general,  will  involve  both  modes  so  that 
each  heave  and  pitch  motion  will  contain  all  the  oscillatory  components  present.  However,  if 
the  coupling  between  the  heave  and  pitch  motions  is  weak,  i.e.,  if  the  contributions  of  the 
A(  .  B  .  and  C'(J  (i  *  j)  terms  in  Equations  (I )  and  <2)  have  small  effect  on  the  coefficients  of 
the  characteristic  Equation  (ft),  then  useful  insights  can  be  obtained  by  both  examining  the 
uncoupled  heave  and  pitch  characteristics  and  comparing  them  with  the  corresponding  coupled 
values.  Furthermore,  for  the  present  case,  even  when  the  coupling  was  very  important,  it  was 
found,  for  SWATH’s  4A  and  4B.  that  the  period  of  one  of  the  modes  was  always  close  in 
value  to  the  uncoupled  heave  period  and  was  easily  identifiable.  It  was.  therefore,  meaningful 
in  this  report  to  refer  to  “heave”  and  "pitch"  modes.  The  uncoupled  periods  are  readily 
shown  to  be 


for  heave 


for  pitch 


(II. I) 


(11.2) 
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These  have  meaning  only  when  the  quantities  under  the  radicals  are  positive,  i.e.,  the  mode 
is  oscillatory.  This  was  always  true  for  tne  heave  equation  and  almost  always  true  for  the 
pitch  part. 

By  analogy  with  the  simple  oscillator,  each  mode  can  be  thought  of  as  having  a  damping 
ratio  fn  which  is  defined  in  terms  of  the  pair  of  roots  describing  that  mode.  Thus,  if  X,  and 
X2  (Equations  (8))  are  identified  as  the  roots  of  the  heave  mode  then 


X,  +  x2 

’  Vvl 

and  similarly  for  the  pitch  mode 


(12.1) 


(12.2) 


Another  measure  of  the  degree  of  stability  is  the  time  that  must  elapse  for  the  disturb¬ 
ance  motion  of  each  mode  to  double  or  halve  itself.  In  the  case  of  oscillatory  modes  it  is  the 
envelope  of  the  disturbance  that  doubles  or  halves.  It  is  easy  to  show  that  this  time  is  given 
by 

T«/2i  or  ■!<*»«  Bn  2/|XnR  i  (12.3) 

Clearly  as  a  stable  mode  becomes  progressively  less  stable,  the  value  of  T*1^'  becomes  larger 
and  larger.  It  is  generally  the  root  with  the  smallest  real  part  that  is  critical. 

The  corresponding  values  of  the  damping  ratio  and  the  time  to  half  amplitude  for  the 
uncoupled  heave  and  pitch  modes,  respectively,  can,  of  course,  be  written  much  more  simply 
in  terms  of  the  coefficients  of  the  equations  of  motion.  They  are  given  as  follows  for  ready 
reference. 


B33 

fhO*  —  l<M»AJ)>C„r’'1  (13.1) 

fpO  =  —  Ul5  +  A5J)  Cjjl  1/2  (13.2) 


II 


tU/2*  = 
'hO 


2(M  +  Aj3  )  fin  2/Bjj  for  fh0  <  1 
2(M  +  A3j  )  fin  2  j 


Bjj  -  {  Bjj  -  4(M  +  A33)C 


„l »] 


(13.3) 


for  fh0  >  I 


2(IS  +  A55 )  fin  2/B55  for  rp0  <  I 


n'o2'  ■ 


2(1,  +  A,,  )  fin 


■/ 


B 


'55 


B^j  4(1,  +  A55 )  (  ,5 


1/2 


(13.4) 


for  fPo  >  1 


In  addition  to  the  transient  response,  the  frequency  response  of  the  craft  to  sinusoidal 
waves,  which  is  often  represented  by  response  amplitude  operator  (RAO),  is.  ol  course,  also 
very  much  affected  by  the  stability  roots.  Much  insight  can  be  obtained  from  an  examination 
of  the  equations  for  RAO  expressed  in  terms  of  the  stability  roots,  especially  at  the  important 
encounter  frequencies  near  resonance. 

If  we  replace  the  right-hand  sides  of  Equations  (I  )  and  (2)  by  the  wave-exciting  terms 
j:  t,iu)i  anj  [.^uwi  respectively,  then  we  can  show  that  the  heave  and  pitch  motion  are. 

respectively,  given  by 


U  - 


I  j  e,w' 


a  (iw  X,  )  (iw  X;)(iw  X,  >  (iw  X4  ) 


(14  1  ) 


*5 


_ _ 

a  (iw  X(  )  (iw  X, )  (iw  X, )  (iw 


(14.2) 


where  a  is  defined  by  Equation  (7.1 ),  and 


Fs(  w2AJ5  +  iwBj,  +  ^  35 ) 

Fj(  w2  A5j  +  iwBjj  +  CSJ) 

and  X(.  i  =  I, ...4  are  the  stability  roots. 

As  will  be  seen  later,  the  heave  mode  roots  X, .  X,  for  SWATH's  4 A  and  4B  are  complex 
conjugates  for  all  the  conditions  of  interest.  Furthermore,  this  mode  is  significantly  less 
stable  than  the  pitch  mode.  We.  therefore,  take  a  closer  look  at  Equations  (14.1 )  and  (14.2) 


Fj  =  F,  |  w*(I5  +  Aj, )  +  iwBJ5  +  f  55  | 
Fj  =  F5  j  w2(M  +  A,3)+  iwB33  +C„  j 
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at  wave-encounter  frequencies  near  the  damped  natural  frequency  of  the  heave  mode.  We 
rewrite  these  equations  as  follows 


u 


2  -  \  x 

OU  AiA? 


I  2 


=  square  of  undamped  natural  frequency  ot  heave  mode 


(X.  +  X,  )* 

►2  = 

3  4X.X2 


(X3  +  X4  )2  \ 

^XjX,,  / 


square  of  damping  ratio  of  heave  (or  pitch)  mode 
(Hquations  (12  1  )  and  ( I  2.2)) 


t  (cj)  and  e2(w)  =  arguments  of  heave*  and  pitch-mode  terms,  respectively,  in 

1  the  denominator 


It  can  be  seen  from  the  previous  equations  that  the  term  in  the  first  bracket  in  denominator 
D  can  have  a  sharp  minimum,  depending  on  the  magnitude  ot  the  damping  ratio  f  3 .  If  the 
second  bracketed  term  in  D  varies  relatively  slowly  with  frequency,  then  heave  and  pitch 
RAO  might  be  expected  to  exhibit  a  peak  at  the  heave  resonant  frequency  with  character¬ 
istics  mainly  determined  by  the  heave  roots.  As  may  be  inferred  from  Equation  (15).  sharp¬ 
ness  of  the  peak  would  decrease  with  increasing  The  frequency  at  which  this  peak  occurs 
is  given  by 


WH  =  woh  V  1  “£3  '  Xll 

where  X„  is  the  damped  natural  frequency  of  the  heave  mode,  and  wH  is  essentially  the 
resonant  frequency  for  the  heave  mode.  It  is  clear  from  the  previous  equation  that  the 


1  2*f 
.»  -tL 


1/2 


for  £,  < 


vr 
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resonance  frequency  is  smaller  than  both  the  undamped  and  damped  natural  frequences. 
Substituting  this  expression  for  u>„  in  the  bracketed  terms  of  Equations  (15.1 )  and  (15  -) 
leads  to  the  following  equations  for  the  approximate  peak  response. 


F', 


U  =  T 


J  l(UM  t,(CJ)  €j(w)) 


a  Du 


(15.3) 


*5  = 


FS  i(u>t  -  fjtwo 


(15.4) 


where 


D„  =  ^on  ^  "*? 


OH 


XJ  X4 


(1  2f>) 


4-  4 


^5h 

X,X4 


(1 


-ri > til 


1/2 


(15.5) 


and 


< 


yr 


-is  also  an  approximate  indication  of  the  magnification  factor  of 


The  quantity  - 

«3  /I  t  j 

,|,c  response  relative  .. the  static  response  <«-0l  which,  according  to  Equations  1 1 5 1.  would  e 

1 

approximately  proportional  to 

aw0H  XJX4 

As  has  been  noted,  the  foregoing  statements  are  reasonably  correct,  provided  the 
numerator  and  the  terms  contain, n,  X,.  X,  in  the  denominator  are  no.  overly  sensitive  to 
frequency  in  the  vicinity  of  heave  resonance.  This  condition  is  satisfied  in  par.  by  the  vehicles 

considered  here,  especially  for  the  condition,  in  which  the  pitch  mode  ***°™"” 
stabk  than  the  heave  mode.  It  should  also  be  noted  from  Equations  .14)  that  pitch 
motion  will  not  have  as  sharp  a  peak  as  the  heave  motion  since  the  pitch-exc, tattoo  term 
tends  to  approach  a  minimum  in  the  vicinity  of  the  heave  resonance  .hen  the  heave  damptng 

is  low  see  equation  for  Fj  in  Equation  414.:.  Similar  considerations  to  those  given  previous  y 

for  the  heave  mode  would  also  apply  to  the  pitch  mode.  Thus,  if  both  the  pitch  and  heave 
modes  have  low  damping  ratios,  then  the  frequency  response  .ill  show,  m  general,  two 
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resonant  peaks,  corresponding  to  resonant  encounter  frequencies  in  heave  and  pitch  modes. 
This  is  true,  provided  the  resonant  frequencies  are  reasonably  well  separated.  When  they  are 
close  together,  and  the  damping  is  low.  a  much  smaller  minimum  in  the  denominator  of 
Equations  (15*  can  occur.  Although  one  may.  therefore,  expect  a  sharper  resonance  peak, 
this  can  be  modified  appreciably  by  the  fact  that  F'  and  F',  also  may  get  smaller  near 
resonance. 

One  motion  criterion  of  great  importance  is  represented  by  the  magnitude  of  the  rela¬ 
tive  bow  motion  with  respect  to  the  wave  surface.  This  can  be  obtained  by 

R  =  it,  M,  tji 


iiu  “V  t  ♦  uii 

where  77  =  A  e  is  the  wave  elevation  at  the  bow.  and  is  the  x  coordinate 

of  the  bow 

We  can  express  the  relative  bow  motion  in  a  nondimensional  form  as 


R  !  ^Ms  HU)**  ♦  U!l  1 

X  1  ~a~  “  c 


(lb) 


The  value  of  R/A  depends  on  the  amplitudes  and  the  phases  of  heave,  pitch,  and  wave 
motions.  However,  if  we  choose  fins  which  provide  good  stability  and  small  heave  and  pitch 
amplitudes,  we  are  more  likely  to  have  smaller  relative  bow  motions.  From  Equations  (15.3) 
and  ( I  5.4)  we  have  seen  that  the  amplitudes  of  heave  and  pitch  motions  1$,  I  and  If,  I  at  the 
resonance  of  the  heave  mode  increase  rapidly  as  the  damping  ratio  decreases.  The  same  is. 
oi  course,  true  for  the  pitch  mode.  By  increasing  the  size  of  the  stabilizing  aft  fins,  it  will 
be  seen  that  adequate  pitch-mode  damping  can  be  readily  obtained  while  the  same  is  not 
necessarily  true  for  the  heave-mode  damping.  Thus,  our  aim  should  be  directed  toward 
obtaining  a  fin  configuration  which  would  also  provide  suitable  heave  damping  and  a  natural 
frequency  so  that  the  heave  amplitude  at  the  heave  resonance  becomes  minimum. 


NUMERICAL  RESULTS 

Input  data  necessary  to  obtain  the  stability  roots  Xn  from  Equation  (b)  are  given  in 
Table  2.  The  subscripted  coefficients  B*  are  obtained  from  available  experimental  data  for 
submarine  models,  and  the  remainder  of  the  subscripted  coefficients  A  and  B  are  obtained 
from  potential  flow  theory.5  including  free  surface. 


Pien,  P.  (  .  and  < .  M  U>e. 
Symposium,  Pans  (1172). 


‘Motion  and  Rcsutanro  of  a 


Low- Walcrpljnc- Area  Calamaran."  The  Ninth  Naval  Hydrodynamic 
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TABLfc  2  INPUT  DATA 


NondimensionaliMtion 

Factor 


SWATH  4A 


SWATH  4B 


P 

g 

V  (Displaced  Volume) 

w 

L  (Overall  Length) 
A_ 


M 


>*p 


(GM)V 


'5 

A33 

A35  A53* 


'55 


33  f  b33<‘>d 
•'L 

!  =  J<-K)  b33 

■'l 

1/ 


<°! 

J35 


101 


(x)dx 


B56=|  */b33,1,|dx 


BH 


b;6 


b;5 


B53 


'LCK 


V/L, 

V 

Li 

pv  L3 

pV 

PV  L, 
PV  L3 

pvv^h 

pVv^ 

pvliN/5^ 


*L3 

2 

—  L3 
2 

—  L4 
2 

P  3 

—  L3 
2 

1 

1 

1 


2.0 

32.17 

83,027  ft3 
208.3  ft 

239.7  ft 
4  732 

-0.0495 

0.133 

00625 
0629 
0  00517 

0  0585 

0.111 

0  000582 

000567 

0.00494 

00010 

00004 

0  006 
3.0 

-85  0  ft 
80.25  ft 


2  0 
32  17 

83,027  ft3 
247.3  ft 

275.9  ft 
5  746 

-C  0475 
0.236 

00625 

0.668 

0.00409 

0.0587 

0.0885 

0000369 

000486 

0.0036 

0.00078 

0.00026 

-0.004 

3.0 

-101  .C  ft 
101.9  ft 
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One  of  the  conditions  for  stability,  according  to  Routh  in  Equation  (9).  is  that 

e  >  0 

On  substituting  Equation  (7  5)  for  e.  we  obtain 

1171 

From  Equations  (3)  and  (17.1 )  we  have  a  speed  condition  for  stability  for  the  bare  hull  case 
as  follows 


which  states  that  a  sufficient  condition  for  instability  is  that  the  speed  be  greater  than  the 

(\5<v, 

right  hand  side.  Since  AC.Mj  >  — - for  SWATlEs  4A  and  4B.  as  can  be  shown  by  sub- 

stituting  the  values  from  Table  2  for  the  coefficients.  Equation  (17.2)  may  be  simplified  to 


Tit  is  amounts  to  approximating  the  inequality  (17.1  »  by 

c55  >  0  <17  4' 

The  computed  values  of  U„  for  SWATII’s  4 A  and  4B  are.  respectively.  25.1  and  3b. 0 
knots.  The  larger  value  of  U„  for  SWATH  4B  is  mainly  attributable  to  the  larger  C,MS.  which 
is  58.4  feet,  compared  to  27.7  feet  for  SWATH  4A.  Values  of  U()  found  in  this  manner  are 
very  close  to  those  obtained  from  the  coupled  equations. 

To  maintain  a  reasonable  stability  to  40  knots,  both  ships  would  need  a  stabilizing  fin  or 
Pms.  To  determine  the  appropriate  fin  size,  a  first  approximation  for  the  lower  bound  is  made 
from  Equation  (3  12)  to  satisfy  Cjj  >  0  at  40  knots.  That  means 
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N 

A  iJMj  -  U2  -y-  l’m;  4  U‘c,.,8,  clm  >  0 

fl 

For  N  *  2.  i.e.,  one  fin  for  each  hull  of  identical  size,  and  6,  *  S2  <  0,  we  have  for  the 
lower  bound 


CIS1  ^LOl 


4u2L3M;  -  A  GMg 
p  U2  IS,  I 


(18) 


for  U  =  40  knots.  Although  the  fin  size  obtained  in  this  manner  would  be  inadequate,  it 
can  be  used  to  represent  the  lower  bound  of  the  ranee  of  fin  sizes  to  be  investigated. 

With  the  aid  of  earlier  estimates  for  the  SWATH  4  fins,  obtained  in  this  manner,  estimates 
tor  SWATH  4 A  and  SWATH  4B  showed  that  the  following  fins  are  reasonably  good  sizes  to 
insure  pitch  stability 


Chord 

n 

Location 

SWATH  4A 

11.9 

14.3 

Inboard  side  of  each 
hull  at  0.84  L 

SWATH  4B 

10.0 

12.0 

Inboard  side  of  each 
hull  at  0.85  L 

It  should  be  stressed  that  what  is  really  being  optimized  in  this  investigation  is  the 
product  c, s,  CLai.  Clearly  this  quantity  may  be  realized  by  an  infinite  number  of  fin  plan 
forms.  The  plan  form  described  previously  is  felt  to  be  conservative  on  the  basis  of  strength 
considerations.  Optimization  of  the  plan  form  for  minimum  drag  for  various  values  of 
cisi  CL<*i  is  bey°nd  the  scope  of  this  investigation;  it  should  be  done  in  the  preliminary  design 
stages.  These  fins  are  designated  1 .0  for  SWATH  4A  and  SWATH  4B,  respectively.  Other 
fins  with  the  same  aspect  ratio  but  different  sizes  are  designated  with  numbers  indicating  the 
projected  fin-area  ratio  to  the  1.0  fin.  For  instance,  the  0.0  fin  means  a  bare  hull. 

For  each  fin.  the  four  roots  of  Equation  (6)  are  found  for  various  speeds,  and  the 
natural  periods,  the  half-decay  time  and  the  damping  ratios  are  computed  from  Equations 
(10)  and  (12).  The  values  are  shown  in  Table  3  for  SWATH  4A  and  in  Table  4  for  SWATH 
4B  for  various  fin  sizes  at  the  speeds  of  30,  35,  and  40  knots. 
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TABU  .<  TRANSIENT  CMARACTI  RISTICS  ()l  SWATH  4A 


Fin  S'/* 

— 

Mode 

Damping  Ratio 

Natural  Period 

sec 

Halt  Decay  Time 

sec 

Aft 

Fore 

35 

40 

30 

35 

40 

30 

35 

40 

knots 

knots 

knots 

knots 

knots 

knots 

knots 

knots 

0 

0 

Heave 

BB 

wm 

9  36  ‘ 

if! 

61  93' 

Pitch 

kb 

KlB 

16  31* 

ECS 

74  96* 

0  85 

mm 

0  70 

0  2? 

9  16 

909 

9  02 

4  88 

4  34 

kb 

080 

0% 

73  51 

31  80 

175  73 

7  60 

7  78 

1  0 

0 

Heav* 

0  75 

0  ?9 

966 

9  78 

4  1? 

7  78 

Ptch 

KB 

0  70 

0  79 

19  05 

71  08 

7  35 

2  07 

i  o  y» 

mm 

0  73 

0?7 

Ci  31 

991 

10  14 

10  44 

4  57 

3  53 

iH 

060 

068 

1/81 

_ 

18  84 

7018 

7  6? 

1  98 

1  10 

0  75 

0  79 

0  3? 

1061 

11  15 

4  46 

3  59 

Pdrh 

0  58 

066 

0  71 

16  63 

16  89 

1706 

7  51 

1  85 

1  15 

0 

Heave 

0  75 

RH 

0  31 

10  68 

Pith 

058 

0  71 

15  5? 

1  TO 

0 

Heave 

0  75 

■ 

mam 

11  01 

11  70 

mm 

imeis 

P'ICh 

0  59 

mm. 

mam 

14  51 

14  11 

■Sa 

WBm 

1  50 

0 

D9 

i  ■ 

1?  14 

1758 

11  40 

15  14 

KB 

17  47 

12  7? 

1  19 

1  01 

20 

0 

Heave 

1?  55 

1?  94 

49  76 

117  88 

Unstable 

P  »Ch 

HI 

1?  74 

11  94 

099 

_ 

0  83 

1  8 

0? 

Heave 

0  70 

0  17 

| 

17  83 

13  51 

1 

,/'V;  ■ 

Pi(Cl> 

0  80 

0  0b 

Kb 

1381 

11  71 

mm 

BIB 

1  6 

04 

0  44 

050 

11  46 

1?  59 

3b 

7  80 

m 

096 

1  0 

76  05 
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Figure  2  shows  the  locus  of  the  stability  roots*  obtained  from  Equation  (6)  for  SWATH 
4A  at  30  knots.  The  blank  circles  with  the  numbers  indicating  the  fin  size  are  for  a  pair  of 
aft  fins.  The  solid  circles  are  the  stability  roots  for  two  pairs  of  fins;  for  each  hull,  one  pair 
is  located  near  the  bow.  and  the  other  is  located  near  the  stern.  Figures  3  and  4,  respectively, 
are  similar  plots  to  Figure  2  for  35  and  40  knots.  As  noted  earlier,  the  vehicle  is  stable  as 
long  as  the  real  part  of  the  root  is  negative. 

The  bases  of  selecting  the  most  desirable  fin  size  are.  as  stated  earlier,  to  insure  stability, 
to  obtain  adequate  damping  for  heave  and  pitch,  and  to  maintain  reasonable  natural  periods. 

As  can  be  observed  in  Figures  2  through  4.  except  for  the  bare  hull  case,  the  aft  fins  ranging 
from  0.7  to  2.0  appear  to  provide  necessary  stability  in  the  pitch  mode.  The  larger  the  fin. 
the  larger  the  absolute  value  of  the  real  root  for  pitch.  This  trend  is  not  true  for  the  heave 
mode.  At  approximately  the  1.1  fin.  IXR I  takes  a  maximum  value  and  decreases  as  the  fin 
size  increases  or  decreases.  From  Equation  (12).  it  can  easily  be  shown  that,  when  the  roots 
are  complex,  the  damping  ratio  is  obtained  by  the  cosine  of  the  angle  between  the  negative 
real  axis,  and  the  line  connecting  the  root  to  the  origin:  vis.,  f  =  cos  7  in  Figure  2  for  the 
damping  ratio  of  the  heave  mode  using  the  1.15  fin;  see  also  Figure  10.  We  find  that  in  the 
case  of  aft  fin  only,  the  1.1  fin  provides  the  maximum  damping  ratio  for  SWATH  4A  between 
30  and  40  knots.  In  this  speed  range,  the  natural  period  for  heave  is  about  10.5  seconds  and 
for  pitch  about  16.8  seconds.  These  periods  correspond  to  encountering  head  sinusoidal 
waves  of  1400  and  2900  feet,  respectively,  at  30  knots.  This  indicates  that  the  natural  periods 
are  reasonably  large,  and.  furthermore,  they  are  well  separated.  The  natural  periods  for  heave 

and  pitch  at  zero  speed  are  9.4  and  15.3  seconds,  respectively. 

The  case  of  using  stabilizing  fins  both  fore  and  aft  was  investigated  to  a  limited  extent. 

As  can  be  noted  from  Figures  2  through  4  for  the  case  of  aft  fins  alone,  increase  of  aft  fin 
size  after  a  certain  size  causes  a  decrease  in  heave  stability  as  well  as  heave  damping  ratio. 

This  fact  implies  that  a  possible  way  of  obtaining  a  large  heave  damping  ratio  for  a  given 
projected  fin  area  is  to  divide  the  fin  area  into  fore  and  aft  fins.  However,  since  the  fore  fins 
generate  a  destabilizing  pitch  moment,  the  aft  fins  should  be  sufficiently  large  to  compensate 

the  destabilizing  moment  induced  by  the  fore  fins. 

The  total  fin  area  was  kept  twice  that  of  the  1.0  fin  but  was  divided  into  fore  and  aft 
fins  of  the  same  aspect  ratio.  The  forward  fin  was  assumed  to  be  located  at  approximately 
the  0.15  L  position  at  the  inboard  side  of  each  hull  and  the  aft  fin  at  the  0.84  L  position  at 
the  inboard  side  of  each  hull.  The  numbers  beside  the  solid  circles  indicate  the  aft  fin  size. 
Hence,  the  size  of  the  fore  fin  is  two  minus  the  number  indicating  the  aft  fin  size.  As  can  be 


•Only  the  positive  imaginary  pail  of  the  complex  root*  are  shown. 


21 


REAL  ROOT 


Figure  2  -  Stability  Roots  for  Various  Fin  Sizes  for  SWATH  4A  at  30  Knots 


Figure  3  -  Stability  Roots  for  Various  Fin  Sizes  for  SWATH  4A  at  35  Knots 
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seen  from  these  figures,  the  heave  damping  is  markedly  increased  by  the  use  of  the  fore-and- 
aft  fin  combination.  However,  the  pitch  stability  begins  to  decrease  with  the  aft  fins  less 
than  the  l  .6  fin  as  speed  increases. 

Since  the  use  of  fore  and  aft  fins  may  lead  to  significantly  higher  drag,  feasibility  should 
be  determined  in  terms  of  the  advantages  they  may  offer  as  far  as  improvement  in  ship 
motions  and  control  is  concerned.  This  is  beyond  the  scope  of  the  present  investigation. 

The  stability  roots  for  SWATH  4B  at  30.  35.  and  40  knots  are  shown  in  Figures  5 
through  7.  For  the  aft  fin  only  cases,  the  size  of  the  aft  fins  investigated  range  from  0.36  to 
2.0  fin.  For  the  fore  and  aft  fin  cases,  the  sizes  of  the  fore  and  aft  fins  are  chosen  like 
SWATH  4A  by  keeping  the  total  projected  fin  area  to  twice  that  of  the  1.0  fin  of  SWATH 
4B  The  numbers  beside  the  solid  circles  indicate  the  size  of  the  aft  fin  of  the  fore-and-aft 
fin  combination.  The  aft  fins  are  located  at  the  0.85  L  position,  and  the  fore  fins  are 

located  at  the  0.12  L  position.  .  .... 

SWATH  4B  without  stabilizing  fins  has  a  speed  ol  inception  for  t  e  pile  mo  e  i 
„V  of  lb  knots.  As  .  an  be  seen  from  Figures  5  through  7.  the  0.0  fm  shows  negative  real 
roots  for  pilch  at  35  knots,  implying  that  the  boat  ,s  stable  in  pilch  even  at  this  high  speed, 
tins  ,s  due  to  the  large  longitudinal  metaeentnc  height  CM,  that  contributes  directly  to 
restoring  moment  of  pitch.  However,  an  increase  in  «M,  decreases  the  natural  penod  of 
push  at  zero  speed.  As  the  speed  increases,  the  destabilizing  Munk  moment  increases  rapidly 
a,  the  square  of  the  speed,  and.  hence,  the  effective  restoring  moment  for  pitch  decreases. 

Tins  results  ,n  an  increase  of  the  natural  period.  Will'  a  stabilizing  fin.  the  loss  of  restoring 
momenl  due  to  the  Munk  moment  is  compensated  to  some  extent,  and  the  natural  period 
does  no,  change  with  speed  too  much.  In  fact  „  ,s  possible  with  a  large  enough  f,n  to  cause 

a  reduction  in  the  natural  period  of  pitch  with  speed:  see  Table  4. 

In  the  case  of  aft  fin  alone,  the  fin  which  provides  a  maximum  heave  mode  damping  is 
found  from  Figures  5  through  7  to  be  the  13  fm  fo,  SWATH  4B.  The  natural  periods  for 
heave  and  P„ch  in  the  30-  to  40-kno,  range  are.  respectively.  0.8  and  14  seconds.  Although 
these  are  reasonably  large  natural  periods,  the  gap  between  the  two  periods  ,s  not  as  large  as 
that  for  SWATH  4A.  The  natural  periods  at  zero  speed  are  0.4  seconds  for  heave  and  - 
seconds  for  pitch.  Tims,  with  the  1.3  Fm.  SWATH  4B  can  widen  the  gap  between  the 
natural  periods  of  heave  and  pitch  as  speed  increases.  If  we  place  mom  emphasis  on  the  gap 
between  the  two  natural  periods  than  on  the  heave  damping,  probably  the  1 .0  Tin  would  e  a 
beer  choice  As  can  he  seen  from  Table  4  for  the  1 .3  and  1 .0  Tins,  the  damping  ratio  tor 
heave  a,  30  knots  reduces  from  O  IK  for  the  1.3  fin  to  0.15  for  the  1.0  fin.  The  correspond¬ 
ing  time  to  half-amplitude  has  also  gone  up  approximately  30  percent,  while  the  di  Itrcncc 
between  the  two  natural  periods  increases  from  4.08  to  6.04  seconds. 
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The  fore  and  aft  fins  for  SWATH  4B  also  demonstrate  an  increase  in  pitch  stability 
and  heave  damping.  The  increase  appears  not  so  pronounced  as  in  the  case  of  SWATH  4A. 
The  combination  of  the  0.4  fin  forward  and  1.6  fin  aft  shows  a  maximum  heave  mode  damp¬ 
ing  ratio. 

In  the  range  from  30  to  40  knots,  the  stabilizing  fins  have  more  effect  on  SWATH  4A 
than  on  SWATH  4B.  This  phenomenon  may  be  expected,  since  in  this  speed  range,  SWATH 
4A  can  maintain  its  stability  only  by  using  stabilizing  fins,  whereas  SWATH  4B  is  stable  to 
36  knots  without  fins  as  a  result  of  the  larger  stiffness  contributed  by  the  lengthened  struts. 

In  other  words,  the  stiffer  the  ship  the  harder  it  is  to  influence  or  control  its  initial  charac¬ 
teristics.  If  the  fins  are  controllable,  the  effectiveness  of  the  fins  may  be  expected  to  be 
significantly  greater  for  SWATH  4 A  than  SWATH  4B. 

The  effect  of  speed  on  the  transient  characteristics,  damping  ratio,  natural  period,  and 
half-decay  time  for  SWATH  4A  with  the  l.l  fin  and  SWATH  4A  with  1.2  fin  are  shown  in 
Figures  8  and  9.  It  is  seen  that  by  installing  the  stabilizing  fins,  not  only  is  the  craft  stabil¬ 
ized  but  stability  is  also  improved  with  increasing  speed.  It  appears  that  SWATH  4A  has 
somewhat  better  overall  stability  characteristics.  Although  this  will  probably  result  in  better 
response  to  the  seaways  of  interest,  it  would  be  necessary  to  determine  this  separately. 

To  study  the  effect  of  fins  on  the  motion  characteristics  of  ships  in  waves,  the  heave  and 
pitch  motions  in  regular  head  waves  were  examined.  Since  our  major  interest  lay  in  deter¬ 
mining  the  critical  motion  behavior,  the  study  of  the  motion  was  limited  to  the  heave  reso¬ 
nant  frequency  region.  A  ship  speed  of  35  knots  was  chosen  lor  this  study. 

If  we  exar.  ine  the  heave  and  pitch  motion  near  the  heave  resonance  given  by  Equations 
(I  5).  we  can  infer  that  to  have  minimum  values  of  Ifcjl  and  l{5l.  we  should  have  minimum 
values  of  IFjl  and  IFjl  as  well  as  a  maximum  value  of  D„  Since  the  effect  of  fins  on  the 
wave-exciting  terms  Fj  and  Fj  is  very  small,  the  minimization  of  Inland  1 15 1  should  be 
accomplished  by  maximizing  the  denominator  DH  on  which  the  fins  have  a  direct  influence. 

We  can  write  DH  given  by  Equation  (15.5)  as 

=  “  W0H  \/*  ^UH  ~  *  <‘U,H  *4  ^ 

where 

UH  =  W0H  >/*  “£}  ^°r  ^ 3  <  r^~ 

V  *- 

is  the  frequency  at  which  the  heave  motion  would  peak,  if  the  magnitude  of  the  term  in 
brackets  were  to  vary  slowly  with  frequency  near  this  frequency.  The  undamped  heave 
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Figure  8  Transient  Characteristics  versus  Speed  for  SWATH  4A  with  l.l  Fin 
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natural  frequency  w0H  and  the  heave  damping  ratio  f,  are  functions  of  the  stability  roots 
X,  and  X2  as  shown  earlier.  The  maximization  of  DH  may  be  pursued  in  terms  of  the  stabil¬ 
ity  roots  of  the  heave  mode  alone,  if  it  yields  a  correct  answer.  For  this  purpose,  we  define 
a  new  denominator,  which  ignores  the  effect  of  the  pitch  mode. 

^ho  =  woh  f 3  \/»  "  fj  =  ^ir  (19) 

In  Table  5  the  values  of  and  D^0  are  shown  for  different  configurations  of  fin  for 
SWATH  4A  and  SWATH  4B.  Fin  configurations  which  are  judged  to  provide  insufficient 
pitch  mode  stability  are  not  considered.  As  can  be  observed  in  Table  5.  the  fins  which  pro¬ 
vide  a  maximum  value  of  DH0  do  not  necessarily  provide  a  maximum  value  of  D„.  This 
implies  that  the  effect  of  pitch  mode  roots  on  the  peak  value  of  heave  amplitude  is  not 
negligible.  As  tar  as  the  value  of  Dj10  is  concerned,  the  fins  that  provide  a  maximum  heave 
damping  ratio  yield  a  maximum  value  of  D|J0.  From  Equation  (15.5)  we  note  that  \}\4  is 
also  an  important  factor  to  be  considered  in  the  maximization  of  DH .  For  those  fin  config¬ 
urations  considered  in  Table  5.  X3  and  X4  are  a  pair  of  complex  conjugate  roots.  Thus. 

XjX4  =  l\/.  and  we  can  denote  this  by  cojjp.  which  is  the  square  of  the  undamped  natural 
frequency  of  the  pitch  mode.  To  best  visualize  the  relation  between  the  stability  roots  and 
the  various  quantities  introduced  in  the  present  analysis,  a  vector  diagram  is  given  in  Figure 
10.  As  can  be  seen,  the  product  of  the  magnitudes  of  the  vectors  from  the  point  iw  to  each 
ol  the  tour  roots  is  simply  DH ;  see  Equations  (14).  Figure  10  thus  shows  the  effect  of  the 
root  distribution  on  the  magnitude  of  DH  . 

It  would  be  still  premature  to  expect  that  the  peak  value  of  the  heave  amplitude  could 
be  a  minimum,  when  the  value  of  DH  is  a  maximum.  We  have  to  include  the  effect  of  the 
wave-exciting  terms  to  judge  which  fin  configuration  yields  the  most  desirable  motion  in 
waves. 

Figure  1 1  shows  the  heave  amplitude  and  the  relative  motion  amplitude  of  the  hull  at 
the  0.06  L  position  for  SWATH  4A  and  at  the  0.05  L  position  for  SWATH  4B  versus 
encountering  period  in  seconds.  These  results  were  obtained  from  an  existing  computer 
program,  which  is  based  on  the  theory  described  in  References  5  and  6.  Four  different  fins 
were  chosen  for  the  motion  study.  These  were,  respectively,  the  best  aft  fins  judged  from 
the  transient  characteristics  the  1.1  fin  for  SWATH  4A  and  the  1.0  fin  for  SWATH  4B-a  fin 
greater  than  the  best  fin.  a  fin  smaller  than  the  best  fin,  and  the  best  combination  of  the  fore 
and  aft  fins. 


Hidlet,  1.  B.  el  al "Ocean  Catamaran  Seakeeping  Design  Baaed  Upon  the  Experience  of  USNS  HAYES,' 
Society  of  Naval  Architects  and  Marine  Engineers,  Vol.  82  (1974). 
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Figure  10  -  Vector  Diagram  in  Stability  Root  Plane 
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As  cun  be  seen  from  Figure  1 1 ,  the  periods  at  which  the  peak  values  of  heave  amplitude 
occur  are  very  close  to  the  heave  natural  periods  given  in  Tables  3  and  4  as  determined  from 
the  stability  roots.  The  periods  at  which  the  peak  values  of  the  relative  motion  near  the  bow 
occur  are  slightly  less  than  the  corresponding  periods  for  heave  motion.  This  is  understand¬ 
able  in  view  of  liquation  (lb).  The  fins  with  greater  heave  damping  ratio  show  a  lesser  degree 
of  rise  and  fall  of  the  motion  amplitudes  near  the  resonance,  which  may  be  expected.  It  is 
interesting  to  note  that  the  configurations  with  the  1.5  and  2.0  fins,  which  were  the  largest 
shown,  had  the  greatest  motion  at  resonance.  This  was.  of  course,  due  to  the  fact  that  the 
heave  damping  ratios,  as  determined  from  the  stability  analysis,  were  only  0.12  and  0.07. 
respectively.  The  values  were  less  than  half  the  value  determined  for  the  best  configurations; 
see  Tables  3  and  4. 

The  effect  of  the  wave-exciting  term  on  the  peak  heave  motion  is  apparent  for  the  0.67 
fin  for  SWATH  4B.  From  Table  5.  the  fin  which  gives  a  maximum  value  for  D„  can  be  found 
to  be  the  1.0  fin.  whereas  the  0.67  fin  shows  a  slightly  better  response  than  the  1.0  fin  as  can 
be  seen  in  Figure  1 1.  However,  in  general,  the  frequency  response  of  the  ships  with  different 
fin  configurations  near  the  heave  mode  resonance  is  relatively  well  reflected  by  the  values  of 
D(J  in  Table  5.  Although  it  is  seen  that  for  SWATH  4A.  the  maximum  value  of  D„  is  obtained 
for  the  0.85  fin.  it  should  be  noted  that  the  pitch  mode  is  beginning  to  become  less  stable 
than  for  the  larger  fins.  It  can  be  seen  from  Figure  4  that  the  configuration  with  the  0.7  fin 
is  on  the  verge  of  instability  in  pitch  at  40  knots.  In  the  selection  of  an  optimum  fin.  on  the 
basis  of  the  present  analysis,  where  many  of  the  coefficients  have  been  estimated,  it  would 
appear  that  a  conservative  approach  should  be  used,  and  that  a  fin  size  yielding  somewhat 
larger  stability  should  be  selected. 

Within  the  context  of  the  present  analysis  and  under  the  consideration  of  pitch  stability, 
motion  characteristics  in  head  waves,  and  the  natural  periods,  it  appears  that  for  aft  fin  alone, 
the  1.0  to  1.2  fins  for  SWATH  4A.  and  the  0.7  to  1.2  fins  for  SWATH  4B  are  the  best 
choices  of  those  considered.  For  the  fore  and  aft  fin  combination,  the  0.4  fin  forward  and 
the  1  .6  fin  aft  for  both  SWATH  4A  and  4B  appear  to  be  the  best  choices. 

CONCLUDING  REMARKS 

Suitable  sizes  for  stationary  stabilizing  fins  for  SWATH’s  4A  and  4B  were  sought.  The 
investigation  revealed  the  following  results. 

1.  Inception  speed  for  instability  is  well  predicted  by  F.quation  (17.3)  derived  from  one 
of  Routh’s  necessary  conditions  for  stability. 

2.  Within  the  context  of  the  present  work  it  appears  that  suitable  aft  fins  would  be  as 
follows: 
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Number 
of  Fins 

Location 

Chord 

ft 

Span 

ft 

Aspect 

Ratio 

SWATH  4A 

One  for 
each  hull 

0.84  L  at  the 
inboard  side 
of  each  hull 

12.0-13.0 

14.4-15.6 

1.2 

SWATH  4B 

One  for 
each  hull 

0.8S  L  at  the 
inboard  side 
of  each  hull 

8.5-11.0 

10.2-13.2 

1.2 

3.  For  a  given  hull  with  aft  fins  only,  there  seems  to  exist  a  certain  size  of  fin  which 
provides  a  maximum  heave  damping.  As  the  size  of  fin  is  increased  beyond  a  certain 
point,  the  stability  and  damping  ratio  of  heave  mode  are  decreased. 

4.  Combination  of  fore  and  aft  fins  demonstrates  the  possibility  of  improved  stability 
and  motion  characteristics  over  aft  fins  alone;  however,  resistance  may  be  increased  due  to 

increase  in  appendages.  Further  analysis  is  required  to  determine  whether  this  would  be 
worthwhile. 

5.  The  larger  fins  introduced  pronounced  coupling  effects  between  heave  and  pitch 
mode.  In  fact,  doubling  the  optimum  fin  size  resulted  in  an  unstable  heave  mode. 

—  6.  Although  SWATH  4B  has  a  larger  inception  speed  for  instability  due  to  her  larger 
GMC,  the  overall  transient  characteristics  of  SWATH  4A  with  an  optimum  fin  are  judged 
better  than  those  of  SWATH  4B  with  an  optimum  fin. 
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